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Abstract 

We study the Casimir effect at finite temperature for a massless scalar field in the parallel 
plates geometry in N spatial dimensions, under various combinations of Dirichlet and Neu- 
mann boundary conditions on the plates. We show that in all these cases the entropy, in the 
limit where energy equipartitioning applies, is a geometrical factor whose sign determines 
the sign of the Casimir force. 
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INTRODUCTION 



The vacuum expectation value (VEV) of the Hamiltonian of a free scalar field in a large 
volume V (so that the allowed Fourier modes tend to a continuum) in N spatial dimensions 
at temperature zero is given by the following expression (given, e.g., in [27|) 



Eo= {0\Hfree\0) = 



(27r)^2c^fc 



N 



(1) 



where ul = k"^ + m^, fc^ = ^^fc: 



i=l 



and we use natural units h = c = ks = 1- This is the zero point energy of harmonic 
oscillators integrated over all momentum modes and over all space. The sum clearly diverges, 
but one may define a regularized subtracted hamiltonian H' = Hfree — Eq. This simple shift 
is equivalent to normal ordering H' =: Hfree '■■ In this way {0\H' \ 0) = for a free field. 
Similarly, once the field is constrained by imposing boundary conditions and thus described 
by a different hamiltonian H, its vacuum expectation value Eq = |0) diverges (|0) is 

the ground state of the constrained hamiltonian). 

It is of physical interest to compute the difference Eq — Eq of the two (divergent) vac- 
uum energies. This was first considered by Casimir (for the electromagnetic field) and is 
therefore called the Casimir energy (and in general the subject is called the Casimir effect). 
As written, this difference is not well defined and thus requires regularization. The properly 
regularized finite part of that series is, by definition, the Casimir energy of the constrained 
system. This is carefully considered in J6|, p 24] and 2l|, p 100] . We perform the regulariza- 
tion via the damping function method [6|. We regularize the energies of the bounded and of 
the unbounded system using the same regularizing parameter A. In Eq.(2) we schematically 
designate this by 

Ecasimir = ( Eq — Eq) . (2) 



reg 



At the end of the calculation we take the limit A = 0. 

The Casimir energy depends on the geometry of the constraints. For example, in Casimir's 
original work he considered two infinite parallel plates with separation d. In this case, 
the Casimir energy, Ecasimir, (per unit area of the plates) is a function of the separation and 
it gives rise to a force on the plates, the Casimir force. 



Ecasimir 



dE, 



Casimir 



dd 



and this leads to a measurable effect, the Casimir effect. 

Various methods may be used to calculate Casimir energies 
method 
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path integration method 
mode summation method j9||, fl§l,fl9| 
a, 0,15], [20], [23 



(3) 



e.g., the Green's function 



22 



14 



dimensional regu- 
etc). There are 



larization method 

several reviews and books on the subject, e.g.. 

In the present work we consider the Casimir effect for a free massless scalar field in 
spatial dimensions due to the presence of two parallel hyperplanes at distance rf, at 

finite temperature. Hence, instead of the vacuum expectation values in the definition of 

the Casimir effect at zero temperature above, we have to consider expectation values with 

respect to the thermal equilibrium density matrix of the system. 
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We shall consider two cases. In the first case the scalar field is constrained to vanish on the 
hyperplanes, imitating the presence of ideal conductors (Dirichlet-Dirichlet (DD) boundary 
conditions (b.c.)). In the second case the field is set to zero on one hyperplane while the 
normal derivative of the field is set to zero on the other, thus imitating the presence of a 
perfect conductor and a perfectly permeable material (Dirichlet-Neumann (DN) b.c). 

In this paper we mostly apply the mode summation method. In the zero temperature 
case we apply the Green's function method as well. Our conclusion is that the DD case gives 
rise to attractive forces between the boundaries while the DN case gives rise to repulsive 
forces, a result which holds for both cases for any dimension and any d > and T > 0. The 
case of Neumann-Neumann (NN) b.c. gives rise to the same effect as the DD case, as we 
shall see in the Green's function section. A case of particular interest is to obtain the force 
per unit area (pressure) on the boundaries in the high temperature limit (namely for T such 
that § > 1 where T^= 

A useful quantity, the Casimir entropy, may be defined 0] , 01 , [Hi , [23] ; |25| • Its sign in 
the high temperature limit determines whether the force on the boundary is attractive or 
repulsive (or zero if it vanishes). We obtain the high temperature limit of the Casimir 
entropy with the following results: For > 1, in the DD case the Casimir entropy is 
negative (attractive force); in the DN case it is positive (repulsive force). For = 1 we 
found that the Casimir entropy is zero in this high temperature limit for both DD and DN 
b.c. (the force between the boundary points tends to zero). Calculations for the Casimir 
energy for the DN case in 3 spatial dimensions were performed in |i26i] via zeta function 
regularization method. The Casimir energy in N spatial dimensions for the DD case was 
obtained in [l|, [l9|. We perform such calculation for the DN case in N spatial dimensions 
via the mode summation method and discuss the Casimir entropy in N spatial dimensions 
for both the DD and DN cases. 



MODE SUMMATION AT ZERO TEMPERATURE 

The dependence of the sign of the Casimir force on the boundary conditions is not yet 
fully understood. It was noted ([ll|, 23]) ^^at for simple geometries in the classical limit 



(i.e., at temperatures wherein energy equipartitioning is applicable) the total Casimir force is 
entropic and the Casimir entropy depends solely on the geometry and the type of boundary 
conditions (by the latter we mean various combinations of Dirichlet and Neumann boundary 
conditions in the geometry of two parallel plates). We are interested in calculating the 
Casimir energy and free energy of a massless scalar field in the parallel plate geometry 
of two hyperplanes located at xat = and = d. Consider N dimensional space with 
coordinate vector x = {xi, ...,X]\f) and a cube of edge length L in this space with faces 
given by Xi = — -j, ^ (i = 1, — 1) and xn = {^^) chosen in order to make 

our volume of interest a finite one. Later we will take L — ^> oo, so that only the boundary 
conditions (b.c.) on xn = 0, d will be important (the reason of treating the coordinate xn 
slightly differently is a matter of convenience). 

The hyperplanes at x^ = and xn = d with b.c. on them divide the cube into three re- 
gions — {^^) < xj^ < 0, < xj^ < d, d < xjy < (-^y^) . We say that our system is subjected 
to DD b.c. if on the hyperplanes x^ = 0,d the scalar field is constrained by the conditions 
4> (ti ^)\xj^=o = ^(^i^)]xjv=d ~ ^- Similarly we say that our system is subjected to DN b.c. 
if the scalar field fulfils the equations (p {t, = dx,^ (f) {ti ^)\xN=d ~ ^- ■^'-'^ convenience 
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the boundary conditions on the faces Xi = ^ {i = 1, N — 1) are taken to be periodic, 
which we may write as (t, x)\^,^_l = (t, (z = 1, — 1) . As for the remain- 

ing faces of the cube {xn = we impose Dirichlet b.c. (p {t, x)\^^^_^l^-^ = 

(p {t, x)\^^^L+d = for simphcity. We expect the details of the boundary conditions on the 
faces of the large cube to be unimportant in the limit L — oo. 

Following our discussion prior to Eq.([2]) ,we define the Casimir energy (Casimir free 
energy) of the scalar field as the regularized difference (in the sense previously discussed) of 
the total mean energy (free energy) of the scalar field for the constrained and unconstrained 
cases. Introducing the Fourier modes of the field, each mode at temperature T contributes 
a mean energy 



(where (3 = and free energy 



^fc(^) = + log (1 - exp {-pnwk)) . (5) 
At zero temperature both expressions coincide and reduce to 

E^{0) = FkiO) = ^nuk (6) 

The angular frequency satisfies uj^ = ck. The solutions of the free massless Klein Gordon 
equation = in the whole space, for DD and DN b.c. on the hyperplanes xn = and 
Xn = d, give rise to the following wave vectors in the region < Xf^f < d : 

DD case: = (^ni, ^nAr_i, ^nAf) 1 

DN case: k = m, f ^ (n^v - |)) J (7) 
rij = {0,±1,±2, ....} for j = 1,...,N -1 and hn = {1,2,...}. 

Similar expressions hold in the other regions. The total mean energy (free energy) of a given 
system is just the regularized sum over the energies (free energies) of the individual modes. 
The Casimir energy is defined by 

{d, T)=Iy1 ^J^. (^) + E ^a:. (T) + E^kz (T) - (T) [ (8) 

I, ki k2 fcs ko 



reg 



and the Casimir free energy is: 



{d, T) = {Y F^,, (T) + Y Fk2 (T) + E (T) - ^ (T) \ (9) 

:i k2 kz ko 



reg 



where the vectors /ci, k2^ k^ correspond to each one of the three regions (the order of re- 
gions is from negative xn to positive) of the constrained system, and ko corresponds to the 
unconstrained system. The sums in Eqs. (IHl) and are discrete. But because L is large 
we can replace by integrals each of the sums whose summation index is multiplied by y- 
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or -^3^. Those — 1 components of k become continuous variables of integration and the 
integrations are performed over the whole space. In the continuum limit the Casimir energy 
(Casimir free energy) in both DD and DN cases is the difference of the total energy (free 
energy) of the constrained system in the volume between the planes and the energy (free 
energy) of the free system in the volume between the planes. The damping function we use 
is the exponential and its explicit form is given in Eq.l ITOi) below. The Casimir energy at 
zero temperature due to DD b.c. is given by 

oo „ 

(Y, / d^-'k^k^ + (z^)2exp(-Av/F + (Hi)2)_ 

n=l'' (10) 

oo 

UkN j d^-^k^¥+^e^^{-\^¥+k^)) 







where k"^ = kl + ... + k'^_^. The damping factor A in Eq. fllUp serves to regularize the otherwise 
divergent integrals; only at the end of the calculation do we take A to zero. For the electro- 
magnetic field we can introduce a cutoff damping factor for purely physical reasons because 
a physical conductor becomes transparent at high enough frequencies, namely frequencies 
higher than the plasma frequency of a given material. For a scalar field we may regard it as 
a mathematical device, since no scalar massless fields were observed in nature so far. The 
Casimir energy due to DN b.c. at zero temperature is obtained from (fTOj) by shifting the 
summation index in the integrand n — > n — |. During the calculation rf-independent terms 
arise, but we omit them since such d-independent terms don't influence physical observables 
such as force and pressure These terms may be interpreted as self-energy of the plates, since 
they are proportional to their area L^~^ . We denote shifted Casimir energy and shifted 
Casimir free energy by the same letters as the unshifted one. The Casimir energy at zero 
temperature due to DD b.c. which we obtain is: 

(rf, 0) = --^^r (^) c (iv + 1) (11) 

(47r) 2 d" 

(which agrees with [l|,[l3])- On the other hand, for DN b.c. the Casimir energy is: 

id, 0) = T^w^r (^) (1 - ^)aN + 1) (12) 

(which agrees with [26]). The Casimir pressure for any b.c. at zero temperature is given by 

P.(.,0,.-^(^). (13) 

We observe that Eq. lITT]) implies attraction between the hyperplanes, while Eq. fll2l) implies 
repulsion. 



MODE SUMMATION AT FINITE TEMPERATURE 

Let us now turn to the temperature dependent case. We are interested to obtain the 
expressions for Casimir energy, Casimir free energy and Casimir entropy in the cases of DD 
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and DN b.c. As in the zero temperature case we consider the differences of energy (free 
energy) of the constrained system and energy (free energy) of the free system with mean 
energy (free energy) per mode given respectively by Eq.(jl]) and Eq.([5]). We may decompose 
the total Casimir energy and total Casimir free energy into a sum of two terms: the Casimir 
energy at zero temperature (which we already know) and a remaining temperature dependent 
term. In order to obtain expressions for Ec^^^ {d, T) and Fc^^^ {d, T) and Ec^'^^ {d, T) , 
p{DN) Poisson summation formula [l2| 



exp (2ninx) = 6 (x — n) (14) 

n=— oo n=— oo 

and our slightly modified version which is obtained by x — x + | 

(-1)" exp i2mnx) = S fx - L - (15) 

n=— oo ra=— oo ^ ^ ^ 

We drop divergent but d independent terms which occur since they don't change the pressure 
on the boundaries. After we perform the integrations and set A = we replace an infinite 
series by another infinite series. This new series expresses the result by known functions 
(McDonald functions K^y [lH) and leads to some simplifications of the expressions, such 
as cancellation of E^ {d,0), and makes easier the calculation of the high temperature limit 
{j3 0). The Casimir energy of a massless scalar field in iV > 1 dimensions at finite 
temperature due to DD and DN b.c. that we obtain is: 



oo 

JV 



m,n=l 

oo 
m,n=l 

(d-dependent divergent terms arise during the calculation of the Casimir energy due to DD 
and DN b.c. for = 1 but they cancel each other at the end. An independent short 
calculation based on ([8]) for the = 1 case confirms that f|T6|) . f|T71) are valid for = 1, 
too.) The double series converges for all /3 > and d > 0, and as /3 — > the series converges 
to zero. The conclusion is that the Casimir energy at high temperature tends to zero in 
accordance with the following heuristic argument [ll|. Each mode carries mean energy given 
by Eq.(jl]). At the high temperature limit, namely for T such that ^ ^ 1 where = |^ 
is the geometry dependent temperature scale, we obtain iJ^ (/3 — > 0) = ^ + O (P), i.e., each 
mode carries the same average energy (which equals the temperature by the well-known 
equipartition property of classical physics). The Casimir energy is the difference of sums 
of mean energies of the constrained and unconstrained systems. Moving the hyperplanes 
xn = 0, d adiabatically apart will change the energy levels but at each step there is one to one 
correspondence between the energy levels of the systems. Hence zero Casimir energy at the 
high temperature limit merely reflects the fact that the number of states of the constrained 
and free systems are equal. Generally the Casimir energy due to equipartition in the high 
temperature limit is of the form 

Ec {d, T ^ oo) = {Nconstrained " A^/ree) X T (18) 
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where Nconstrained and Nfree are given by the integrals of the density of the modes Pconstrained 
and pfree, and correspond to the total number of states of the constrained and free systems. 
In our infinite parallel planes geometry the mode densities are given in ^Jj] and Nconstrained 
and Nfree turn out to be both infinite but equal in the sense that the right side of Eg. (ITS]) 
is zero. 

The Casimir free energy of a massless scalar field in > 1 spatial dimensions due to DD 
and DN b.c is: 







« 9 ™ 1 «-L1 / , " — \ fl I 

2^-^d^-^f3^ + ^ ^-^ 2 \ 1^ J 

m,n=l 



/ \ (19) 

/ Anmnd \ 

— jr-f^iY. 

1=1 
r N-i 

oo ^ (20) 

m,n=l 

(the last two formulas give the correct result for the = 1 case, too (see also below); C (1) 
which is divergent is cancelled by an appropriate term in the double series). The Casimir 
pressure on the boundary for any b.c. is given by 

In the DD case it is negative for every temperature and every d (the McDonald functions, 
Kjy, are decreasing functions of their argument). Hence the Casimir force on the hyperplanes 
= 0,(i is attractive while for the DN case the pressure on the boundary is positive for 
any (3 > and d > 0, which implies repulsive force. 

In the case of one spatial dimension (A^ = 1) there are no transverse modes and the wave 
vector k has only one quantized component. Another way to obtain the expressions for the 
A^ = 1 case is to promote A^ to a continuous variable in (fT6!) . (fTTl) . (fT9|) . (!20|) . represent it as 
l + e(0<e-Cl) and take the limit £ — 0. The explicit expressions (which coincide with 
those obtained by direct calculation) are: 

oo 

^) {d,T) = -pJ2 sinh"' (^) (22) 

n=l 

oo 

^) (d,T) = -^Y. (-1)" si^h"' (^) (23) 

n=l 

oo ^ 

^) {d, T) = -Y: I . (24) 

^ ^ exp (^) - 1 



oo 



— , expf^^l-l 
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CASIMIR ENTROPY AND THE HIGH TEMPERATURE LIMIT 



Once given the Casimir energy and the Casimir free energy of a system, we define the 
Casimir entropy 0] , , [HI , [23] , [24] , [ll] by the equation 

E - F 

5c =^^. (26) 

At the high temperature hmit {(3 — > 0) the terms which contain McDonald functions tend to 
zero; thus we obtain that the entropies at the high temperature hmit in DD and DN cases 
for > 1 are: 

{d, oo) = ^g^C (N) = {d, oo) (27) 

-i3Fr\d, oo) 



(Eg. (1271) coincides with the Casimir free energy in the high temperature limit given in |19|). 
We see that the Casimir entropy depends on the geometry (in our case the separation d) and 
on the type of boundary conditions involved. Its sign is not restricted to be non-negative 
since the Casimir entropy is the difference of the entropies of the constrained system and of 
the free system and thus might be negative, as happens for the DN b.c. Recall that in the 
classical limit the Casimir energy tends to zero due to equipartition, and therefore in this 
case the Casimir force depends on the Casimir entropy only. Hence in the classical limit, the 
sign of the entropy determines the sign of the Casimir force. For the two cases of DD and 
DN b.c, we see that the entropies for the two cases differ in sign, thus giving rise to forces 
of opposite signs. Let us observe that = 1 is the only dimension for which the Casimir 
entropy at the high temperature limit does not depend on d and thus it is an irrelevant 
constant and the Casimir force is zero. In the N = 1 case the boundary consists of two 
points, which do not "produce" enough geometry for the Casimir entropy (which gives rise 
to Casimir force). Another interesting observation which is deduced from Eqs. fll9l) . fl2U]) and 
(pTj) is that the Casimir pressure Pc{d,T) {Pc^^\d,T) or Pc^^\d,T)) has a fixed sign for 
fixed d and any T for > 1. 



GREEN'S FUNCTION AT ZERO TEMPERATURE 

Now we turn to the problem of the validity of the mode summation technique. The 
correlation function of a massless scalar field at two points separated by the hyperplane 

xat = 0, on which the scalar field is subjected to Neumann boundary conditions ( -S^ = 
0) is not zero but is given by 

{T(j){x)(f){y)) = A,{x-y) + A,{x-y) ,^9) 
for X]\f > 0, yN < ^ ' 

where y = {y^yi, ...,yN^i, —yN) and Ai{x — y) = (0| (x) (y) |0) is the homogeneous 
solution of DxAi {x — y) = (in three spatial dimensions Ai (x) = —^P.P.-^ Since 
the field operators at points on opposite sides of the Neumann hyperplane are correlated, it is 
not possible to expand the field into independent eigenmodes living on opposite sides of the 
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Neumann hyperplane. Therefore we should check our assumption of expanding the field into 
eigenmodes separately in each one of the regions in the DN case. In contrast, the correlation 
function at points on opposite sides of the Dirichlet hyperplane {(f)\^^^Q = 0) is zero. In 
this case an expansion of the field into independent eigenmodes living on opposite sides of 
the Dirichlet hyperplane is obviously justified. The expression for the Feynman propagator 
(which is the time ordered VEV of the two point function (0| T(f) (x) (y) \0) = iAp {x — y) 
in spatial dimensions) is given by [3] 

{x,y)N,m = 



where (x — y)^ = riap (x" — y") {x^ — y^) with r]af3 the diagonal metric of the Minkowski flat 

space-time with the following components r/oo = —Vn = ••• = —Vnn = 1- Hj)^^ (x) is the 
Hankel function |2i]. The leading terms of Ap {x, y)j^^^ as m ^ is found to be 



^^- ^r (^) / i_ 



JV-l 
2 



Af (x, ?/)^=i^^_o = ^ ln[m2 (x - yY] (32) 

For > 1 it is m independent and by substituting = 3 into ( l3Ti) we obtain the well 
known result For = 1 we see that as m ^ we obtain an infra-red logarithmic 
divergent term (infra-red divergence). Actually we may obtain Eq. (l32|) if we promote A^ to 
a continuous variable (A^ = 1 -|- e (e <C 1)) in Eq. (l3T|l and keep the leading terms for small 
e (up to d-independent divergent terms) 



Af {x,y)j 



1 



'N>l,e^O 4ni (33) 

'-g/- exp(-| log (x - yf) = ^ log(£2 (x - yf) + O (e) 



We see that Eq. fl33|) coincides with Eq. fl32|) if we take e = m. 

The free Feynman propagator in A^ > 1 spatial dimensions is given by Eq. fl31l) . Green's 
functions for DD, DN, NN b.c. in the region between the planes < x^, yN ^ d are given 
by Eq. flMl) . fl35|) . Eq. fl36l) and in the section "Optical Green's function" (below) we show 
how we construct them. (Eq.(!34l) was previously obtained in jsf.) 



QiDD) ^ 

oo 

(Af {x-y + 2ndz) - Ap (x - y + 2ndz)) 

n=— oo 

G'(^^) (x, y) = 

oo 

^ (-1)" {Af {x-y + 2ndz) - Af {x - y + 2ndz)) 

1= — QO 

(x,?/) = 

oo 

^ {Af {x-y + 2ndz) + AF{x-y + 2ndz)) 



(34) 



(35) 



(36) 
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Once the Green's function with the proper b.c. is obtained, we may use it to express the 
VEV of the energy density of the scalar field in the following way 

(37) 



W = ^0' (x)+(v0(x)) 
I lim^^o[<9xo5yo + ■■■ + dxj^dyj^] {T(f) (x) cf) (y)) 



where in the last line we regularized the hamiltonian density operator using the point split- 
ting = + e^. The VEV of the total energy in volume V is: 

{H)y = j (h) dV (38) 

V 

The Casimir energy density is, by definition 

{h), = {h) - (39) 

where (h) is given by Eq. fl371) for the constrained field and {h) ^^^^ is the energy density of 
the free field. The subtraction in (15^ removes the short distance local divergences in flTTI) . 
which cannot depend on the boundary conditions. The Casimir energy density is the VEV of 
the energy density of the field which is constrained by some boundary conditions, measured 
relative to the energy density of a free field. The Casimir energy is given by 

Ec = ij (x) dv\ +{d independent) , (40) 

\v I reg 

where V is the region between the hyperplanes (0 < xtv < c?), and the (maybe infinite) d- 
independent contribution comes from the infinite region outside the hyperplanes. Henceforth 
we shall drop the (i-independent contribution, which does not affect the Casimir force. In our 
three cases the Green's functions in the region Q < xn , y^ < d are infinite series (sum over 
index n) which contain the terms Aj? (s — y + 2ndz) and A^? (x — y + 2ndz) with different 
coefficients, appropriate to the given case. Substitute the expression for the Green's function 
into Eq. fl37|) . If we interchange summation over n and integration with differentiation with 
respect to the coordinates we find the contribution of each term to the Casimir energy. We 
use the coordinate representation of the free propagator A^ given by fl3T]) . fl32|) . The n = 
term which corresponds to the free propagator is cancelled by Eq. fl39|) . Each "even path" 
(in the next section "Optical Green's function", we will explain the terminology "even" and 
"odd" paths) which is represented by the term lS.p {x — y + 2ndz) (n ^ 0) contributes to 
the Casimir energy 

r 1 

while each "odd path" which is represented by the term Ai? (x — y + 2ndz) contributes 

r(-±i) 1 



where In is the length of the path labelled by n {n may be any integer) and is given by 



(42) 



, _ I \2dn\ for "even path" 

" ~ ^ \2xn + 2dn\ for "odd path" ^^'^^ 
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The overall contribution of the "odd" paths to the Casimir energy is d- independent as one 
may verify for all three cases of boundary conditions and therefore unimportant. Con- 
sequently, summing over the index n we obtain the following well known expressions for 
Casimir energy (for any A^) [l9|, j26| which coincide with Eg. fill p and Eq. fll2p . 

i?f ^) {d, 0) = ^^4^ C (A^ + 1) (44) 

ErHd,0) = -(l-^^Ef^Hd,0) (45) 
£;f^)(rf,0) = Ef^)((i,0) (46) 

OPTICAL GREEN'S FUNCTION 

In this section we are mostly inspired by [3], [3] (which construct Green's functions for 
various geometries via classical trajectories of the rays of light) and Let us consider a 
massless scalar field which is subjected to Dirichlet b.c. on the first hyperplane at x^v = 
and Neumann b.c. on the second hyperplane at a; at = d. Our goal in this section is to show 
how one may write the exact Green's function in the region between the two hyperplanes 
by making classical, geometric optical considerations. We shall express the Green's function 
between the hyperplanes in terms of propagators connecting two points via paths which hit 
the hyperplanes and are reflected (according to the laws of geometric optics). For each time 
that the path hits the Dirichlet hyperplane we multiply the propagator by (—1) while each 
time it hits the Neumann hyperplane we multiply it by 1. Formally we characterize each 

path by an index r = ^o"^'', cr^^^ where and cr^^ are the reflection numbers, the number 
of times a given path r is reflected from the first and second hyperplane, respectively. The 
total number of reflections a^'' + for path r is denoted by |r| . Obviously for every 
classical path which we consider, the following relation holds 



(r) (r) 



< 1. (47) 



It is convenient to separate the paths into "even paths" and "odd paths". "Even paths" 
means paths which are reflected an even number of times from the hyperplanes. For these 

(r) (r) 

paths the sum of reflection numbers, |r| , is even and a}^ = a]^ (The last equahty is 
consistent with (H71) ). "Odd paths" means paths reflected an odd number of times from 
the hyperplanes. For odd paths, by definition, |r| is odd and we may divide them into two 

(r) (r) 

classes according to a}y — a)^ = ±1. The total energy is obtained by summation over all 
the contributions of all paths. The total energy due to "even paths" for which |r| = 2n is: 

Ei^en^ = % [ dVAp ix-y) + 



CO „ „ (48) 

2tJ2 (-1)" / dVAp {x-y + 2ndz) = i En=-oo (-1)" / dVAp {x - y + 2ndz) 

n=l / / 
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The first factor in Eq. (H8|) corresponds to tlie patfi |r| = 0. Tlie factor 2 is due to degeneracy 
of the even paths for which \r\ ^ 0. Together they combine to the r.h.s of Eq. lHHIl . The odd 
paths which first hit the first hyperplane ( the class cr^'' — cr^^ = 1) contribute 



i (-1)""^^ / dVAp {x-y + 2ndz) (49) 
while those which hit the second hyperplane first (the class a^j^ — o"^'' = — 1) contribute 

oo „ 

i (-1)"+^ / dVAp {x-y- 2ndz) (50) 
n=i I 

Together they give rise to 

oo „ 

= ^ E (-1)"^' / ^^^^ - ^ + 2^^^) (51) 

n=— oo y 

The contribution to the Casimir energy of the odd paths is actually zero {E^^^J^^ = 0), as we 
already checked in the previous section, and only even paths contribute to it. It is given by 



oo , 

E(^^)=2 Y {-^T j dVAF{x-y + 2ndz) 



(52) 



V 



(where the prime means taking the sum without the n = term). 

One can also easily obtain Ei^^\ Ec^^"^ which coincide with the known results. 



SUMMARY 

To summarize, we presented in this work calculations for the Casimir energy, free energy 
and entropy of a scalar massless field at finite temperature in the cases of DD and DN 
boundary conditions. We used the technique of mode summation for any temperature 
and also the Green's function method for zero temperature. In the case of DD boundary 
conditions we used the usual Poisson summation formula while in the case of DN boundary 
conditions we had to modify it and wrote it in a slightly different form. Later we used an 
identity for infinite series to express our results by means of modified Bessel functions of the 
second kind (McDonald functions). The new results are the Casimir energy, Casimir free 
energy and Casimir entropy due to DN b.c. At the high temperature limit (namely when 
T ^ Tc) we obtain that Kirchhoff 's law holds, namely, that the Casimir energy tends to zero 
in both cases, thus it is independent of the boundary conditions. Another result which we 
obtain is that at the high temperature limit the force between the hyperplanes is attractive 
for the DD boundary conditions and repulsive in the case of DN boundary conditions. The 
fact that the Casimir energy at high temperature limit is zero but the Casimir force does 
not vanish indicates that the Casimir force in the high temperature limit is purely entropic. 

The results we obtained in the zero temperature by using Green's functions coincide with 
those performed by the mode summation. Further research is needed to understand why 
does the mode summation, which inherently assumes non-correlation between the modes, 
work in the DN case. 
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